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Abstract . An escpressioti is derived for the time transformation t - t, where 
t is coordinate time in the solar system barycentric space -time frame of 
reference and T is proper time obtained from a fixed atomic clock on earth. 
This transformation is suitable for use in the computation of high-precision 
earth-based range and. doppier observables of a spacecraft or celestial body 
located anywhere in the solar system; it can also be used in obtaining com- 
puted values of Very Long Baseline Interferometry data types. The formula- 
tion for computing range and doppier observables, which is an explicit fune^ 
tion of the transformation t- is described briefly. 


vi 
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1 . Introdtaetiom 


This paper derives an expression for the time transformation t- r, 
where t is coordinate time in the aolar system harycentrie space-time frame 
of reference and r is proper time obtained from a fixed atomic clock on the 
surface of the earth. ^ The expression is obtained using general relativity; 
however, to the accnraey of the retained terms, it is consistent with all 
viable relativistic theories of gravitation. The expression for t - r is suitable 
for use in the computation of high-precision earth-based range and doppler 
observables of a spacecraft or celestial body located anywhere in the solar 
system, It can also be used in obtaining computed values of Very Long Base- 
line Interferometry data types. The expression for t - t can he used in orbit 
determination programs in which the motion of bodies and Light is repre- 
sented in the solar system harycentrie space -time frame of reference with 
coordinate time t as an independent variable. The errors in computed range 
and doppler observables due to neglected terms in the expression for t - t will 
not exceed approximately 0.62 m and 2.4 X 10 ^ m/s, respectively, per 
astronomical unit of range to the spacecraft. These figures apply specifically 
for two-way tracking data (transmitted and received at the same tracking 
station on earth). 

An expression for t' - r, where t* is coordinate time in the heliocentvie 
space-time frame of reference, was previously obtained by Moyer (1971). 
However, a term affecting three-way tracking data (transimitted at one station 
on earth and received at a second station) was inadvertently omitted. This 
previous expression for t' ^ t was obtained by a straightforward integration of 
the differential equation for dt/dt*. Thomas (1975) has shown that the use of 
integration by parts and a first-order expansion of the gravitational potential 
simplifies the derivation and provides a clearer understanding of the physical 
origins of the various terms. The present derivation for t - t, where t is 
coordinate time in the solar system harycentrie space-time frame of refer- 
ence, uses the method of Thomas (1975) and produces an expression that 
includes all of the terms previously obtained by Moyer (1971), with minor 


Coordinate time t is proportional to proper time that would he read hy an 
atomic clock at infinite distance from the solar system and at rest relative 
to the solar system barycenter. 
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changes in the coefficients, and the above-mentioned three-way term. The 
coefficients obtained in this paper are slightly more accurate than those pre- 
viously obtained. Th.:; differences are due to different truncation errors for 
the two methods. In addition, six new periodic terms are obtained. Two of 
these are due to long-period variations in the gravitational potential at the 
earth due to Jupiter and Saturn. The remaining four terms are due to the 
offset of the solar system harycenter from the center of the sun* Without 
these four terms, the new expression for t - t applies {with slightly reduced 
accuracy) for coordinate time t^ in the heliocenitrie space-time frame of 
reference . 

The differential equation for dr/dt is developed in Sections 2 and 3. It 
is integrated to give an intermediate expression for t - r in Sections 4 and 5. 
The periodic terms of this intermediate expression, which are integrals or 
dot products of position and velocity vectors, are converted to sums of 
sinusoidal functions in Section 6. These terms are eolleeted in Section 7 to 
give the final expression for t - t. Auxiliary equations for eopiputing the 
arguments of the periodic terms of t - T are also given. The final expression 
for t - T and the auxiliary equations give t - r as a function of time and the 
earth -fixed coordinates of the atomic clock. Section S gives estimates for 
errors in computed range and doppier observables due to terms neglected 
in the final expression for t - r. The effects of the retained terms of t - r on 
these observables are also given. 

Section 9 gives an alternate expression for t - T which is a function of 
the position and velocity vectors of the atomic clock and the major bodies of 
the solar system. This expression is the intermediate expression for t^-r 
given in Section 5 with the three terms which are expressed as integrals 
replaced by functions of position and velocity vectors. In certain circum- 
stances, it may be desirable to compute t - t from position and velocity 
vectors using the equation given in Section 9 instead of computing it as a 
function of time using the formulation of Section 7. 

The notation used is defined in the text. However, the definitions of 
the symbols used globally throughout this paper are repeated in Appendix A, 
Numerical values are given for those parameters which appear in the final 
expression for t - r. Appendix B gives equations for computing range and 
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d.Gppl6r observables. These eguations are explieit finiGtions of the titn.e 
transform^ation t - T . 


2 . Differ eatlal Eqoation for dT/dt 


The invariant interval ds between two events with differences in their 

12 3 4 

space and time Goordinates of dx , dx , dx , and dx is given by 


d'S ' 


g. .dx'’ dx^ 


( 1 ) 


where the repeated indices are summed over the integers 1 through 4. The 

matrix of eoeffieients g^^^ is the metric tensor, obtained by solving Einstein's 

field equations. A solution for the case of n slowly moving bodies in the 

weak field approximation is the n-body metric tensor of Eddington and Clark 
2 

(1938). This solution can be applied to the solar system. For this appliea- 

12 3 

tion, the coordinates x , x , and x of this solution are nonrotating rectan- 
gular components x, y, and z, respectively, of position relative to the 

4 

solar system baryeenter; the coordinate x = ct, where c is the speed of 
light and t is coordinate time. In order to distinguish it from coordinate 
times of other solutions of the field equations, it will be referred to as 
coordinate time in the solar system barycentric space -time frame of ref- 
erence. An interval of proper time dr recorded on an atomic clock is 
related to the interval ds along its world line by 


d T 


ds 

G 


( 2 ) 


Hence, (1) and (2) relate an observed interval of proper time to the changes 
in the space and time Goordinates of the atomie cloek. 

For the purpose of obtaining an expression for t - t, where T is proper 
time recorded on a fixed atomie clock on earth, and t is coordinate time in 
the solar system barycentric frame of reference, the n-body metric tensor 


This is the solution previously obtained by de Sitter (1916) except for a cor- 
rection to one of his terms. The n-body metric tensor corresponds to the 
n-body Lagrangian given by Eq. 3.3.37 of Infeld and Plebanski (i960). The 
n-body metric tensor (with, reversed sign conventioh) may also be found in 
Eq. 39.63 of Misner, Thorne, and Wheeler (1973). 
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is sia'bstituited. into (!) and terms are retained to order (1/c)*^, givdng 

ds ~ ( ^ ^ c dt - (dx + d 7 + dz ) (3) 

where U is the Newtonian gravitational potential at the atomic clock, com- 
pnted nsing the positive sign convention (i.e., U = - d )• Let the velocity of 
the atomic clock relative to the solar system baryeenter be denoted by s, 
which is defined by 


•2 

s 



(4) 


Subs titia ting (4) and (2) into (3) gives 



Expanding and retaining terms to order 1/c gives 


(5) 


dr 


dt 



(■ 6 ) 


4 3 

The neglected terms of t ^ r due to the neglected 1/e terms of (6) have a 

- 12 

maximum magnitude of about 10 s. They affect computed range and dop- 
pler observables by a maximum of about 10 ^m and 10 ^^m/s, respectively. 
Equation (3) is the "Newtonian” approximatioin to the n-body metric {see 

4 

paragraph 39*7 of Misner, Thorne, and Wheeler, 1973). All viable 
reiativistic theories of gravitation have the same metric to this level of 


^JDue to terms neglev^-ted in (3) and in expanding (5). 

Equation (3) implies the spaceiike sections of simultaneity t = constant to be 
flat Euclidean spaces to first order; the coordinates x, y, and z io them, 
can properly be thought of as Euclidean coordinates, and all of the results 
of classical Euclidean geometry may be used, e. g. , triangle theorems and 
vector addition laws for slowly moving objects. 


4 
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approximation. Menee, the expression for t-r obtained from (6) is consistent 
with all viable relativistic theories of gravitation. 

Equation (6) is the basic differential equation relating observed atomic 
time T to coordinate time t in the solar system barycentrie frame of refer- 
ence. The rate dr/dt of an atomic clock relative to uniform coordinate time 
is a function of the Newtonian gravitational potential U at the clock and the 
solar system barycentrie velocity k of the clock. 

Equation (6) will be used to obtain an expression for coordinate time t 
minus proper time obtained from a fixed atomic clock on earth. In the 
expression to be obtained, the proper time will be specifically International 
Atomic Time TAI disseminated by the Bureau International de I'Heure (BIM). 
The TAI second is the SI second {International Syst'^m of Units), which is the 
duration of 9 192 631 770 periods of the radiation eorresponding to the transi- 
tion between the two hyperfine levels of the ground state of the cesium- 133 

atom (Meeh'tly, 1969). This is the observationally determined average length 

3 

of the second of coordinate time t. Referring to Equation (6), periodic 

. 2 

variations of U and s from their average values produce periodic terms in 

• 2 

t- T. However, the average values of U and s for a fixed atomic clock on 
earth are positive, and the average length of an interval of proper time dir 
is less than the corresponding interval of coordinate time d,t. Thus, the 
average length of the second of proper time r is greater than th'e length ©f 
the second of coordinate time t- Proper time t in {6) does not correspond 
to TAI because the average length of the TAI second is equal to the length 
of the second of coordinate time t. Let TAI be denoted by T*, which differs 
from T in the length of the atomic second. A modified form of (6) will be 
obtained which applies for atomic time r^-. The average length of the r* 
second must equal the length of the second of coordinate time t. Hence, the 


Interpolation of the lunar ephemeris with an observed longitudie of the moon 
gives an "observed" value of the LHdependf^.^U!t variable, coordinate time t, 

The average number of cycles obtained from a cesium atomic clock per 
second of coordinate time t was obtained by counting cycles of a cesium 
atomic clock between two observations of the moon and dividing the observed 
number of cycles by the diffefenee of the two "observed" values ei coordi- 
nate time t. 
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differeatiai equation relating and t must satisfy the condition that the 
average value o£ dT^ equals dt. 

In (6), dr is obtained as 


dr = 


dN 

n 


( 7 ) 


where dN is an observed number of cycles obtained from an atomic clock 
and n (in units of cycles per second) is a conversion factor from cycles to 
seconds of atomic time. The value of n corresponding to (6) is that value 
which results in the length of the second of proper time equalling the length 
of the second of coordinate time when the atomic clock is an infinite distance 
from the solar system and fixed relative to the solar system baryeenter. 
Equation (6) can be rewritten as 


dr _ _i_sf_ U-U 1 _ 

^ ~ 2“ Z 2 " 2 " 2 2 

e e c e 


( 8 ) 


6 . ' 2 

where U is the time average value of U at the atomic clock and s is the 

time average value of s for the clock, Ignoring 1/c"^ terms, this may be 

written as 


dt 



dr 


R. 1 
2 " 2 



- U I S ^ - 3 ^ 

2 2 2 
c e 


(9) 


Substitute (7) into (9), and let atomic time be obtained as 


dr 



( 10 ) 


The time avera.ge value of U is 


U = lim 

X— *co 



U dt 


6 
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where the Gonversion factor n’!' is given by 


Then (9) is given by 



XJ J_ 

2 “ 2 2 / 
c c / 


( 11 ) 


I 


ill = 1 u-n J_ 

dt 2 " 2 2 

e e 


( 12 ) 


I 


The average value of dr=!= equals dt and hence atomie time is the mathe- 
matieal representation of International Atomic Time TAI. 

• 2 

Perio<iic variations of U and s from their average values result in 
periodic variations in atomie time t-'f- relative to coordinate time t in the 
solar system barycentrie frame of reference. Some of the periodic terms of 
t - T* for atomie clocks at various locations on earth are not in phase. Thus, 
one is tempted to conclude that it is not possible to synchronize a world- wide 
network of atomic clocks which read atomic time 7"i=. However, Section 3 
shows that such a system of atomic clocks can be synehronized with light 
signals . 

In the conversion factor n'^ given by Equation (11), the average value 
U varies with the location of the atomie clock on earth. However, the only- 
term of U which varies significantly is the gy vitational potential U(E) due 
to the earth. Similarly, the average value s varies with the location of the 
atomic clock on eart^ only because of the variation in the square of the geo- 
centric velocity v due to the earth's rotation (v = uoj, where u is the distance 
of the atomie clock from the earth's spin axis and uj is the inertial rotation 

rate of the earth). The only significant variable part of n=!= is thus 

' 2l 2 7 

-n U(E) -t (1 /2) V /e , The quantity U'{E) + (1 /2)v , the sum of the earth's 

gravitational and centrifugal potentials, is constant on the geoid (mean sea 

level). The conversion factor n* thus varies with altitude above the geoid; 

it increases at the approximate rate of 1.00 x 10 ' cycles per second per 

kilometer. If the defined length of the TAI second (9 192 631 770 cycles 

from a cesium atomic clock) is taken to apply at mean sea level, the 

conversion factor n^ (cycles from a cesium atomic clock per second of TAI) 
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9 192 631 770 + 0.001 00 k 


(13) 


is 


given 


where h is altitude above the geoid in kilotineter s . 

TAI is obtained in praetiee as a weighted average of times obtained from 
atomic clocks located at various altitudes. All of the clocks have the same 
conversion factor from cycles to seconds of atomic time. Prior to comput- 
ing the weighted average, time obtained from each clock is corrected in 
value and rate (Bureau International de I'Meure, 1975). The rate corrections 
remove the variations in the clock rates due to the differing aititudes above 
mean sea level and other causes. Hence, in the determination of TAI, the 
altitude -dependent conversion factor n* is used implicitly, but not explicitly. 

In (12), the average lengths of the TAI and coordinate time seconds are 
the same and are equal to the "observed" length of the coordinate time second. 
Hence, by using (12), the estimated length of the second of coordinate time 
becomes its adopted length. This does not produce any errors in observed 
minus computed residuals for tracking data obtained from a spacecraft, if 
the n-body epfaemeris (for the planets, sun, and moon) is fit to observations 
using the expression for t - T'f' obtained from Equation (12). Also, the observed 
minus computed residuals for the planets, sun, and moon are not degraded. 


3 . Synchronization of TAI Atomic Clocks 

This section shows that a world- wide network of atomic clocks which 
use the conversion factor n=i' given by (11) and read atomic time = TAI 
can be synchronized by light signals. This conclusion is reached by con- 
sidering the following expression for the interval ds which applies in a local 
region of the geocentrie inertial (i.e., n 0 nr 0 tatin.g) frame of reference: 


ds 


2 


2ii.(:e) ^ 

2 

c 


c dt r- 


2.^2 ^ 2 , 
-t dy -t dz ) 


(14) 


where U(C) is the gravitational potential due to the earth, t is coordinate time 
in the geocentric inertial frame of referenee, arid x, y, and z are rectangular 
components of position in this frame. Equation (14) is a modified form of (3) 


h^bodxjcibmy of the 
PAGE IS POOR 
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applied to a local region of the geoeentrie inertial frame of reference. The 
modification consists of replacing the total gravitational potential by the 
potential due to the earth only. From the Frineiple of Equivalence, the 
accelerating earth produces a gravitational field which cancels (in a small 
region surrounding the earth) the gravitational field due to the sun, moon, and 
other bodies which aeeelerate the earth. Thus, in the local geoeentrie line 
element (14), the gravitational potential due to all bodies other than the earth 
does not appear. 

2 

Substituting (.2) into (14) and retaining terms to order 1/c gives 


dT _ , U(E) 1 

It - ^ ■ ~T~ ‘ 2 le j 

c 


(15) 


where s is the velocity of the atomic eloelt relative to the geocentric inertial 

frame of reference. The sum of the gravitational and centrifugal potentials, 

■ ? 

U(E)-t(l/2)s , IS constant on the geoid and thus the rate of proper time t varies 
with altitude above the geoid. The interval dr is obtained as 


dr 


m 

m 


( 16 ) 


The value of the conversion factor m (cycles per second of atomic, time) which 
corresponds to Equation (15) is that value which results in the length of the 
second of proper time T e quailing the length of the second of coordinate time 
t for a fixed atomic clock at infinite distance from the earth. If the length of 
the atomic second is changed so that 




m 

ra^ 


(17) 


where 


m^' 



e 


1 > 

z-r> 


e 


(18) 
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equation {iS) rediaees to 



U'Sittg the altitude -depeneient Gotiversion factoar results in atomic 

time T5i=, which runs at the rate of coordinate time t in the geocentric inertial 

frame of reference. The constant m in (18) can he selected so that m* = n" 

given hy (11) at mean sea level. However, the altitude -dependent terms in 

(18) and (11) are proportional to m and n, respectively, which differ by terms 
2 -3 

of order 1/c . Thus the factor 1.00 X 10 in (13) will agree with the corre- 
sponding factor obtained from (18) to about 7 digits. The maximum fractional 

-21 

difference between m^ and n^ is about 10 , which is negligible in relation 

to the current stabilities of atomic cloclcs. Thus, to sufficient accuracy, 
atomic time r* in (19) is the same timie scale as t * in (12). 

From (19), atomic cloehs at varying altitudes which read r-l' have the 
same length for the second of atomic time (in seconds of coordinate time fc). 
These clocks could be synchronized with a master clock by using iight signals 
transmitted by way of an equatorial synchronous satellite. The Newtonian 
light time (calculated in the geocentric inertial frame of reference) from the 
master station to any other station is constant. The reLativistie eontrihution 
to the light time is due to the mass of the earth only; it is less than 0. 2 ns and 
constant. Hence, intervals of reception at the various stations will he 
identical to intervals of transmission at the master station. Since the various 
clocks have the same length for the atomic second, the transmission intervals 
can be synchronized with the seconds pulses of the master clock, and the 
reception intervals will he synchronized with (hut out of phase with) the seconds 
pulses of the clock at the receiving station. In practice, only one signal is 
required to synchronize each clock with the master clock. Upon reception of 
the signal, the clock at the receiving station is set equal to the known trans- 
mission time plus the light time calculated In the geocentric inertial frame of 
reference. The signal can be transmitted by way of any spacecraft in the 
earth's vicinity or hy lunar bounce. For other methods of synchronizing 
atomic clocks on earth, see Thomas (1975). 
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4. Integration of Differetiitlal Equation 


Eqiuatlon (12) relates an infinite simal interval of atomic time r* 
obtained from a fixed atomie clock on earth to the corresponding interval of 
coordinate time t in the solar system harycentric frame of reference, the 
Newtonian gravitational potential U at the clock, its time average value U, 
the square of the solar system harycentric velocity of the clock i , and its 
time average value s^. 

This section integrates Equation (12) to obtain an interme'diate 
expression for t- The terms of this expression, which are integrals or 

dot products of position and velocity vectors, are converted to sums of 
sinusoidal functions in Section 6, 

This section gives the magnitudes of several discarded terms of t - r 
and their maximum effects on computed range and doppler observables. 

The effects that are range-dependent were evaluated at a range of 50 astro- 
nomical units. The most significant of these errors are included in the error 
summary in Seetion‘8, 

In the derivation of the expression for t - t the effects of solid earth 
tides, polar motion, and nutation are ignored. The maximum effects of 
these phenomena on two-way or three-way range and doppler observables, 
due to changes In t - are given in Table I: 

TABLE I 

Effects on range and. doppler observables due to changes in t - 


F henome non 

Bange, m 

Doppler, m./s 

Solid earth tides 

10'^ 

lO"”^ 

Folar motion 

le"^ 

lO""^ 


-1 

-4 

Nutation 

10 

10 ^ 


In the evaluation of the terms of t - t--' in Section 6, the effects of precession 
are accounted for. In the final expression for t - r=i=, the coefficients of some 
of the terms are explicit functions of the earth-fixed coordinates of the 
atomic clock. Continental drift will result in a slow drift in the values of 
the s e c o o r dina te s . 
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In (12), U - U exelndes the eonstanit potential dale to the earth and can 
be expressed as 


U-li = (U^-U^)p +(U 


E P 


where and are the gravitational potentials due to all bodies except 
the earth, evaluated at the loeation of the atomic clock and at the center of 
the earth, r espeetively. The subscript p indiGates that only the periodic 
terms of the quantity are to be retained. 


In the following, r|, and r| will refer to the position, velocity, 

and acGeieration vectors, respectively, of point i relative to point j. The 

quantity s| is the velocity of point i relative to poii j. The components of 

rj are space coordinates in the solar system bar'y ^entric spaee^time frame 

of reference. The dots denote differentiation with respect to coordinate 

. 2 

time t in this frame of reference. In (12), s can be expressed as 


.2 

s 



( 21 ) 


where A, E, and C refer to the locations of the atomic clock, the center of 
the earth, and the solar syst'=“m barycenter, respectively. Then 



( 22 ) 


Substituting (20) and (22) into (12), d'enoting r =l= by t, and integrating 
from an Initial epoch (subs crip t o) to the current time t gives 


t “ r = (t - t) 

o 








e 



TT a. * 

- “e 



dt 


p 


(23) 


12 
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In tke first integral of ( 23 ), is nearly constant and periodic 
variations in tbe square of this quantity can be ignored. In the second 
integral, the potential difference may be expanded to first order; 



where the gradient (VU)-rj, is evaluated at the center of the earth. The 

neglected higher-order terms in this expansion produce terms of t- r with 

magnitudes of 10 or less. They affect computed range and doppler 

"•3 - S 

observables by less than 10 m and 10 m/s, respectively. In the second 
integral of (23), the dot product of the two velocity vectors is integrated by 
parts : 


I (^E ■ ^a) p * = "7 (iS -J^a) j ^ JJ^A ■ 4)p '*■' 


Since terms of order 1/c arc ignored, the acceleration of the earth 
can be approximated by its Newtonian value: 




(26) 


and the second term o.: (25) cancels (24). This eaneellation of the potential 
variation earth's acceleration, is in accordance with the 

Principle of Equivalence , According to that principle, in a small region 
surrounding the freely falling earth, the sum of the gravitational potential 
due to the bodies accelerating the earth and the inertial potential due to 
the accelerated earth is constant to a high degree of accuracy. Substituting 
(24) through (26) Into (23) and discarding the term gives 
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t - r = (t - r )q 



( 27 ) 



e 


•C E \ 
) 


t 

t 


Tke derivatioTi of (27), starting from (23), is due to Tliomas (197 3). His 
Equation (13) e or re spends to (27) above. 

In the Integral in (27), the gravitational potential at the earth dne to 
the moon can be ignored. The neglected periodic terms of t - r are smaller 
than 3 ns and affect competed range and doppler observables by less than 
Q. 1 m anid 0. 2 x 10 m/s, respectively. The potential in (27) ean be 
expressed as 



( ’^E - 


+ U 




where U_ and U„ are gravitational potentials at the earth and eartli-m^oon 
barycenter B, respectively, excluding the terms dne to the earth and 
moon. The potential differenee may be ex,panded to first ord'er: 






where the gradient is evalnated at the earth^moon barycenter. The 

neglected terms in this expansion prodnee terms in t- t which are smaller 
- 1 2 

than 10 s and .iffeet eompnted range and doppler observables by less 

_ 4 _ 9 

than 10 m and 10 m/s, respectively. In (27), the sqnare of the solar 
system barycentrie velocity of the earth ean be expressed as 


14 
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(30) 


(4f = (iE+iB)-(iE*i5) 


f (=2)' 


^ 


i® + 

-E 


(•S) 


The periodic terms of j~ can be neglected, Tbey produce terms 

of t - T smaller than 10 ^*^s whick affect computed range and doppler 

- 3 - 8 

observables by less than 10 m and 10 m/s, respectively. The dot prod- 
uct in (30) can be expressed as 





(31) 


Ignoring relativistic terms and assuming the earth and moon are located 
at their baryeenter gives 



(VU)^ 


(32) 


Substituting (28) - (3 2) except the term 


(^ 2 )^ 


of (30) into (27) gives 


t - r = (t - r). 



(33) 
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wkere excludes the potential due to the earth and moon.. The caneel- 
lation of the terms of t-r doae to (29) and the second term of (31) is in 
ac e or d.axiG e with the Principle of Equivalence. This eaneellation is 
associated with the earth-moon barycentric frame of reference and is 
analogous to the previously described cancellation in the geocentric frame 
of reference. 

In (33), the second term of the integral can be expressed as 






(34) 


where S refers to the center of the sun. The last term of the integral is 
negligible. The second term is integrated by parts; 



(35) 


Since terms of order 1/e are ignored, the integral in (35) can be 
evaluated using Newtonian accelerations. Also, the earth and moon can 
be considered to be located at the earth-m'oon baryeenter . The result is 

dt (36) 

P 



WRODUCIBILITY OF TIE] 
' ‘NAL PAG-E IS POOF 
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wtiere the summation over i ineludes eight planets and the earth - moon 

S 

baryeenter. The quantity?^ is the heliocentrie position veetor of planet i, 
r^ is its magnitude, and is the gravitational constant of the planet (see 
Appendix A). When i refers to the earth - moon barycenter, is the suim 
of the gravitational constants of the earth and moon. Assmning that the 
planets move in the ecliptic plane, (36) is given by 



(37) 


where r is the magnitude of , H is the true longitude of the sun rneasured. 

at the earth _ moon barycenter, and is the heliocentric true longitude of 

planet i. Both of these angles are referred to the mean equinox and ecliptic 

of date. When i refers to the earth - moon barycenter, ^tnd 

i-H . ^ ±180*. 

1 


The first term of the integral in (33) can be written as 


7i: 


(u^,) dt = 

B p 


r 

UgCS) 


- 


dit + 




where U_(S) and U_(i) are the contributions to the gravitational potential at 

M ID 

the earth - moon barycenter due to the sun S and a planet i, respectively. 
The sumimation over i ineludes all planets except the earth and excludes 
the moon. For an outer planet L, 
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{H~L) 


, ^ 3 ^ ,2 

1 ( 7 :) (£-1) + €os2(£-^.) 


*^{t) - 1 ( 1 :) 


dit 


( 39 ) 


where all terms to order (r/r^) have heen retained, and the planet i is 
assumed to move in the eeliptie plane. When i refers to an Inner planet, r. 
and r are interchanged in (3 9). 


For an outer planet, the second term of (39) cancels the term of (37) for 
that planet. The first term of (39) is retained far Jupiter and Saturn only. 

The remaining terms of (3 9) are ignored for all outer planets. For an inner 
planet, there is no eaneeilation between (39), after interchanging r and 
and the term of (37) for that planet. All terms of (39) are ignored for 
Mercury and for Venus. The terms of (37) for Mercury, Venus, and the 
earth - moon baryeenter are also Ignored. Substituting the surviving terms 
of (34), (35), (38) and (39) into (33) and repiacing solar system barycentrie 

. Q 

velocity vectors with heliocentric velocity vectors plus r_ g gives 

t - T = (t - f) 

o 




e 




'P 

, 1 

/ , C S \ 

t 

^ 1 

/.C 

E\ 

?: 

c 

(-S ■ -b) 

t 

0 

e 


^a) 


t 

t 


( 40 ) 
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wkere |j.j and are gravitational eonstaMs for Jupiter and Saturn, 


respectively, 
Tlie term 


and r j 


and 


are the heliocentric distances of these bodies. 


M\ 

1 ‘ I 


is negligible and has been omitted from (40). The neglected terms in 
Equations (37) and (39) result in a maximum error in t - r of about 1.3 |j.s; 
the corresponding errors in computed range and doppler observables can 
be up to 1, 1 m and 0. 3 x 10 m/s, respectively. 


For an outer 
at the sun due to the 


the first term of (39) is the gravitational potential 

2 

(divided by c ). The remaining terms are thus 


the potential at the earth - moon barycenter U^(i) minus the potential at 
the sun XJe(i)- The second term of (39) is the dot product of the gradient 
of the potential due to the planet (evaluated at the center of the sun) and 
the heliocentric position vector of the earth - moon barycenter. The 
cancellation of this term with the term of (37) due to the same planet Ls in 
accordance with the Principle of Equivalence. In a small region 
surrounding the sun, the variation in the gravitational potential due to a 
planet Ls cancelled to high accuracy by the inertial potential arising from 
the sun's acceleration due to that planet. However, the cancellation of the 


potential difference U^(i) - Ug(i) to an outer planet i is only 
approximate as evidenced by the highei* -or der terms in (39). For an inner 
planet, the cancellation of this potential difference is not even approximate; 
the term of (37) for an inner planet does not cancel any term of (39) when 


r and r. are interchanged in (39). 
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liiitial Ycdiaeg 


Tile purpose of tkis seetLon is to d'etermiue the value of (t-r)^ in 
Equation (40). The subscript o indicates that this quantity is to be evaluated 
at the initial epoch t^. The value of (t - t )^ is station -dependent; that is, it 
(depends upon the location of the atomic clock on earth. However, it will be 
seen that (t - r)^ minus the initial values of all of the periodic terms of (40) 
has the same value for all stations; it is (denoted by 


The eonelusion that the constant does not vary with the location 

of the atomic clock on earth follows from a consideration of the Lorentz 
transformation. One of the four equations of this transformation is: 



(41) 


In this application, t and t' refer to coordinate time in the solar system 
barycentric and geocentric inertial frames of reference, respectively. 

The cc^ordinate x' is the component of the geocentric position vector of the 
atomic clock along the instantaneous (direction of the earth’s velocity 
vector relative to the solar system barycenter. The magnitude of this 
vector Is s. 

Section 3 showed that atomic clocks on earth which read International 
Atomic Time T-'' (denoted as r starting in Seetion 4) can be synchronized. 

It was also sh(E>wn that this time scale is identical with eoortilinate time in 
the geocentrie inertial frame of reference. Hence, for two sueh 


JFL Technical Memorandum 33-7§b 


synehronized clocks, the interval At' hetween corresponding ticks is zero, 
and (41) gives the interval of coordinate time t between the two "synchronized" 
ticks : 



(42) 


where Ax' is the difference of the x' coordinates of the two atomic clocks. 
Thus, fixed observers in the solar system baryeentric frame of reference 
do not agree that the two atomic clocks on earth are synchronized. The 
difference in the definitions of simiiltaneity in these two frames of reference 
is dne to their relative velocity s. 

The transformation t - r given by (40) mnst be consistent with (42), 

This allows the relation between valiaes of AT^ for different stations to be 
determined. L>et ft - t)^ and (t - t)^ refer to valnes of t - rfor two different 
stations on earth at a given epoch r , Since the atomic clocks at these 
stations (which read r) are synchronized, the difference (t - t)^ - (t - r)^^ is 
the interval of coordinate time t in the solar system baryeentric frame of 
reference between corresponding ticks of the two atomic clocks on earth. 
This interval mnst be eqnal to At given by (42). The sum T of terms 
fonr and eight of (40) excluding the values of the functions at the initial 
epoch is 


T 


2 

c 



4 ) 


(43) 


The value of T for station 2 minus the value for station 1 is equal to the 
approximate form of (42). In order for (t - t)^ ~ (it - r) ^ to be equal to this 
same quantity, the term (t- r)^ of (40) minus the initial values of terms 
four and eight must be the same constant for all stations. The initial 
values of the remaining terms of (40) are not station-dependent; they can 
be subtracted from this constant to give a second constant which is the 
same for all stations. This second constant is AT^ , previously defined 
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as (t - t )^ mirnis the values of all periodie terms of (40) at the initial epoch 

t . Using this notation In (40) and expressing term four as the sum of two 
o 

terms gives 

t-r = AT^ 



where the eonstants of integration for the three ind'efinite integrals are 
zero. 


All terms of (44) except the first term are periodic. Terms four, five, 
and nine represent periodic variations in the rate of a specific clock on earth 
relative to uniform coordinate time t in the solar system harycentrie frame 
of reference. The remaining terms represent periodic variations in the 
rates of all clocks on earth. At a given epoch, the difference in the sums of 
terms four, five, and nine for two different stations on earth represents the 
difference in the definitions of simultaneity for observers on earth and 
observers fixed in the solar system baryeentrie frame of reference. 

The expressions for t - r given in this report apply for ideal International 
Atomic Time r , which in general differs from actual International Atomic 
Time 'r(A). If an e3ipression is desired for t - r{A), terms representing 
T' - T(A) must be added to the expressions for t - r given in this report. 

It was announced at the Ihth General Assembly of the International 
Astronomical Union (held in Grenoble in August 1976) that the second of TAI 


■tWKOi)i:ciEn,n-Y of the 
'’^(.TNAL P&GS IS POT® 
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- L2 

is shorter than the SI second at mean sea level hy 1. 0 ± 0,2 x 10 s. A 

n 

resolutiGn was adopted to increase the length of the TAX second by 
1. 0 X 10“^^s at Jarmary i, 1977, 0^ TAX, Using r = r (A) at this epoch, 
the approximate expression for t - r(A) which applies from January 1, 1972, 
0^ TAX (when the TAX time scale was officially adopted) to January 1, 1977, 
oil tax is: 


r - r(A) =-1x10 ^^t 

where t (a negative number) is seconds past January 1, 1977, 0^^ TAX, The 
a priori estimate for r-f(A) which applies after this epoch is zero. 

The International Astronomical Union has adopted”^ an atomic time 
scale which uses the SI second at mean sea level and has the value TAI 
+ 0. 0003725 days exactly (32, 184 s exactly) at January 1, 1977, 0^1 TAX. It 
has further stated that there be only periodic variations between this time 
scale and coordinate time in the solar-system barycentric space-time frame 
of reference, lienee, for t = r (A) at January 1, 1977 0^ TAX, the Inter- 
national Astronomical Unioin has adiopted 32. 184 s (exactly) for the value of 
the constant 

The value of AT^ does not affect observed minus computed residuals 
for tracking data obtained from a spacecraft if the n-bod'y ephemeris used 
was fit to observations using the same value for AT^. The value of AT^ 
is an adopted constant used in the expressioin for t - t; there is no error in 
AT^ which contributes to errors in eonciputed range and doppler observables. 


7 


See Transactions of the International Astronorruical Union, Vol. XVI B (1977). 
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(n . EvaltaatioH of Periodic Terms 


This section converts the p>eriodie terms in Eq-uation (44) for t - t to 
sTiims of sinusoidal fiinetions . The various as sum'ptiQus made and the terms 
which are discarded are stated but are not justified in this section. However, 
the resulting errors in computed range and doppler observables are given in 
Section § along with the effects of the retained terms of t - r. The final 
expression for t - rand the equations for evaluating the arguments of its 
periodic terms are given in Section 7. 


Several of the periodic terms in (44) contain dot prodiucts of position and 
velocity vectors, AH of the velocity vectors are inertial velocity vectors: 
the inertial velocity vector gives the velocity vector of point i relative to 
a nonrotating coordinate system centered at point j. In evaluating the dot- 
produet terms, the components of each of the two vectors are referred to a 
nonrotating coordinate system which is aligned with the instantaneous position 
of a rotating coordinate system. It should be noted that the velocity vectors 
do not give velocity eomponentc- relative to the rotating coordinate systems. 
The two vectors required for the evaluation of a particular term are referred 
to either the mean earth equator and equinox of date or the mean equinox and 
ecliptic of date. By using these coordinate systems, the final expression for 
t - r aecGunits for the precession of the earth's equator. 

6.1 ANNUA E TEHM 

The second term of Equation (44) is 


(t - r)^ 


c2 



UptS) 4 (4) 



(4S) 


where the constant of integration is zero. This term is evaluated assuming 
that the heliocentric orbit of the earth-moon baryeenter is an ellipse with 
semi-major axis a, eccentricity e, radial coordinate r, and eccentric 
anomaly E. The gravitational potential at the earth^moon baryeenter due to 
the sun is 


U 


1 




(46) 


JPL Technical Memorandum 33-786 



where is the gravitational eonsfcant of the sun. The square of the helioeentrie 
veXocit’/ of the earth -moon, haryeenter is 

(®b) "" ^ 

where and are the gravitational constants of the earth and moon, 
respectively. The inverse of r is given by 

J. ^ l + .leos E (48) 

r a r 

Snibstituting (46), the approximate form of (47), and (48) into the integrand of 
(45) and retaining only the periodic terms as specified by the subscript p 
gives 






2p e 

^ cos E 

r 


(49) 


The derivative of E with respect to t is 


E 



(SO) 


Snbstitnting (49) into (45), multiplying the integrand by d E/dt and dividing it 
by the approxim.ate form of (50) gives 


(t-r)2 



sin E 


(51) 


Substituting numerical values from Appendix A gives, in units of seconds, 


(t - r)^ = 1.658 X 10“^ sin E 
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The eoefficieat is given to four digits heeanse e is eonstant to ahont that many 
digits from 1950 to 2 000. The error dne to this approximation is inelnded in 
the error summary in Section §. The approximation of the factor 
by in (47) and (50) is valid since these two quantities agree to six digits. 

The term (51) was obtained by Moyer (1965). It was also obtained nsing 
an alternate derivation by J, D- Anderson (nnpnblished) in 19^>4. The dei^i’^s-- 
tion of (51) does not involve an expansion in powers of e. Ciemenee and 
SzebeheLy (1967) nsed snch an expansion and obtained a series of terms, the 
first of which is (51) with E replaced by the mean anomaly M. Their com- 
plete series does convert to the form (51), The annnal variation in the rate 
of atomic clocks on earth corresponding to the term (51) of t - t was obtained 
by Aoki (1964). 

6.2 MONTHEY TEEM 

The tbird term of Eqnation (44) is 


(t - T)^ 


2 

e 





which can be written as 


(t - r)^ 


_ I 

C^'(1+|J,) 




where 



is the geecentrie position vector of the moon and 


( 53 ) 


(54) 




(55) 


Both vectors in (54) are evaluated assuming circular orbits, and the 5° inelina- 
feidh of the lunar orbit to the ecliptic is ignored. With rectangular components 
referred to the mean equinox and ecliptic of date, these vectors are given by 
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• B 

-S 


and 


-sin. L 
eos L 
0 


(S6) 


m 

-M 


eos d 
sin d 
0 


J 


M 


(57) 


where 


I 

e 




a 


M 


circular orbit velocity of earth-moon baryeenter relative 
to the sun (see A-p).pendix A) 

earth -moon -barycentrie mean Longitude of the sun, referred 
to the mean equinox and ecliptic of date 

semi-major axis of geocentric orbit of the moon 


(L = geocentric mean longitude of the moon, referred to the 
mean equinox and ecliptic of date 


Substituting (§6) and (57) into (54) gives 


/I 1 ^ 

(t - t) = —Y" sin D 

e (l+|i) 


(S§) 


where D, the mean elongation of the moon from the sun, is given hy 


D = f - L 


(59) 


Substituting numerical values from Appendix A gives 


(t - r)j 


1. 548 X 10"^ 


sin P 


( 60 ) 
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The error sucninary In Seetion 8 gives the errors in cotnipirted range and 
doppler observables dne to truneating the coefficients of the periodic terms 
of t - T . 

6.3 DAII.Y TEEM SUE TO MONTHLY MOTION OF EAETH 
The foiirth term of Eqnation (44) is 


(t - r)^ 



(■ 61 ) 


which can be expressed as 


(t - 


1 

C^(l+p) 




(62) 


The orbit of the moon is assumed to be circular and in the ecliptic plane. The 
vectors in {62), with rectangular eomponents referred to the mean equinox 
and ecliptic of date, are given by 


'-sin ([1 





(63) 

u cos ^) 

- 


u sin (®j^+ cos e 

+ V sin e 

(64) 

■ u sin (©,, \ ) sin 6 

' M 

+ V cos € 



where 

= circular orbit velocity of the moon relative to earth (see 
Appendix A) 

e = mean obliquity of the eeliptic 

“ inclination of the ecliptic plane to mean earth equator of 
date 
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mean sidereal time 
eqninox of date 


Greenwich hour angle of mean 


M 


and n, v, and X are earth-fixed coordinates of the atomic clock. The 
coordinate u is the distance from the ea.rth's spin axis in kilometers, v is 
the distance north of the earth's equatorial plane in kilometers, and X is the 
east longitude. Suhstituting (63) and (64) into (42) and using trigonometric 
identities gives 


(t - ^>4 




n 

2 


(1+eos e ) sin (@j^+X-(X) 


(65) 


(1-cQs e) sin (0j^-t-X + ([) 


-f V sin c- cos C / 


The second and third terms of this expression are neglected. 


From pages 73-74 of the Explanatory gupplemei^ to the Ephemeris 

196 1), the relation hetweeti mean sidereal time , and universal time UT is 

M 


®M UT + Rtj - 12 


( 66 ) 


where UT refers specifically to Uioiversal time UTl, hours past midnight, 
and Ey is the right ascension, measured from the mean equinox of date, of a 
fictitious point on the equator. The adopted expression for H , with units 
changed from hours to degrees, is 


= 279‘’41'27'.'54 + 129 602 768'.' 13 T^ H- 1’1'393 5 


( 67 ) 


where T,, is the numher of Julian centuries of 36 525 days of UTl elapsed 

U lil 

since January 0, 1900, 12' UTl. From page 98 of the above reference, the 
mean longitude L of the sun, referred to the mean equinox and ecliptic of 
date, is given by 
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L = 279“41'4§'.'04 + 129 602 768'.'13 t + 1V089 (68) 

where T is the number of Julian centuries of 36 52 5 days of ephemeris time 
(Goordinate time t) elapsed since January 0, 1900, 12^ ephemeris time. The 
constant term of is 20'.'6 smaller than the corresponding term of L because 
Kjj. is corrected for stellar aberration. The difference between and L 
varies from about 21" at 1950 to about 23" at 2000. Because of this small and 
nearly constant difference, in {'66) is approximated by Li which gives, in 
radians , 

= UTl + L - TT (69) 

If UTl Is given in seconds past January 1, 1950, 0 UTl, the angle UTl in 
radians is computed from 

UTl (radians) = 2 tt j j (70) 

( L J decimal part ) 

The calendar date for the reference epoch for UTl (seconds) in (70) is 
arbitrary, but the time of day must be 0^ UTl. The effects of the approxi- 
mation in (69) are included in the error summary given in Section §. 

Substituting (69) into the first term of (65) and using (59) gives 

s^(H-cos e) 

(t - r). = ~ g — u sin (UTl + \-D) (71) 

2e (1+p) 

Substituting numerical values from Appendix A gives 

(t - r)^ = 1. 33 X 10'^^ u sin (UTl + K-D) (72) 


30 
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6.4 DAILY TI:RMS AND SMALL ANNUAL TERM DUE TO ANNUAL MOTION 
OF EARTH 

TLe fifth term of Equation (44) is 


(t - r)g 


/.S E\ 

N • ^a) 





(73) 


The vectors in (73), with reetaugula-r eomponeuts referred, to the mean earth 
equator and equinox of date, are given by 



-sin (£ -Y) 
cos (i -Y) eos € 
cos (f - Y) sin e 


'S 



u sin 

V 


(74) 


(75) 


whe r e 


velocity of earth-moon baryeenter relative to sun 


i - earth -moon-harycentrie true longitude of the sun, referred 

to the mean equinox and ecliptic of date 

Y = elevation angle of helioeentric velocity vector of earth - 

moon baryeenter, measured from the transverse direction 
(normal to radius) 


From Broucke (1974), 


■ /I i 2 
- S ( 1 - T® ■ 

e 4 


cos M+ve cos 
4 


(76) 


where all terms to order e have been retained, and M is the mean anomaly 
of the helioeentric orbit of the earth -moon baryeenter. The angle (. ~Y 
can be expressed as 


JFL Technical Memorandum 33-786 


31 


i - y 


L + (f-M) - y 


(77) 


where f is the true anornaly ©f the helioeentrlG orbit of the earth-moon bary- 
center. From Herrick (1971), 


V = tan 


-1 


e sin E 


/u?" 


e sin E 


(78) 


where the approximate form is correct to order e . From Smart (1 


sin E 


sin M + sin 2M + 


(79) 


and 


f-M 


5 2 

ie sinM + -^e sin 2M + . 


(80) 


Safes titniting (78) - (80) into (77) gives 

i - y = E + esinM+^e^sin2M (81) 

which inclades all terms to order e . Safes titating (74) and (75) into (73), 
using trigonometric identities, and then sabstituting (69), (81), and (76) 
gives 


(t - r) 


5 


s (1+cQS e )^ 1 2 % 2 

■ — — — 2 (^"J^ 2M) 


X sin 


l+X.-e sin M - sin 2M) 


s ( 1 -cos t )a 

■ ' — (1+e cos M) sin (tJ'T 1 + \ +2 L+e sin M) 
2e 


s (sin e )v 
c 


cos L 
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where all terms to order e , e , and e have been retained in terms one, 
two, and three, respectively. J^xpanding and retaining terms to the same 
order of e gives 

s (1+cos € )n , ~ 

(t - T)g = — (1-ie ) sin (UTl+K) 

2 c ^ 

1 2 

+ e sin (UTl + \-M) - -^e sin (UT1+X+2M) 

O 

+ |-e^ sin (UTl+\ -2M) 


s (l-eos e - 

sin (UT1 + K+2L,) 

2g L 


+ e sin (UT1 + \H-2L+M) 


s (sin e )v 

- — — — 2 — ~ ^ (® 3 ) 


Snbstitntittg nnmerieal valnes from. Appendix A gives 

(t - r)^ = 3. 17679 x iO"^® n sin (UTl + X) 

-12 

+ 5 . 3 12 X 1 0 12 sin (UT1+X.-M) 

_ 14 

- 1. 1 X 10 m sin (UT1 + K+2M) 

+ 1.00 X 10“^^ n sin (UT1+K-2M) 

- 1. 3677 X 10"^^ n sin (UT1+ \+2 l) 

- 2 . 2 9 X 10 n sin (TJT1+ X+2L+M) 

- 1.3184x10 V cos I, (84) 
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The term of (83) and (§4) with argument (UTl+\) was first obtained by- 
Anderson (1968); the expression for the eoeffieient of his term differs slightLy 
from that given in (83). The term with argument UTl + K + ZM is not included 
in the final expression for t-T; the effects ©f deleting this term are included 
in the error summary in Section 8. 

6.5 LONG -PERIOD TERMS DUE TO JUPITER AND SATURN 


The sixth term of Equation (44) is 

where id-,, is the gravitational constant of Jupiter and is the helioeeutric 
J J 

radial coorcainate of Jupiter. The constant of integration is zero. The seventh 
term of (44) is (85), with each subscript J replaced by the subscript SA, 
which refers to Saturn. 


The 
variable . 


inverse of r is given by (48), with a subscript J added to each 

J 

Substituting 1/i’j into (85) and discarding the constant term gives 


" "»6 = “ r 




cos Ej dt 


( 86 ) 


Multiplying the integrand by dEj/dt and dividing it by E ^ given by 


E 


J 


1 




g i ve s 


>* - ’’'4 


sin E J 

c S 


(87) 


( 88 ) 


reproducwity of tee 
-■• 'ITMAL PAGE IS POOR 
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is the circular orbit velocity of Jupiter relative to the sun given by 


wli e r e 


J 



(S9) 


Substituting numerical values from Appendix A into (89) and (88) gives 




5.21 X 10"^ sin E 


J 


The corresponding term for Saturn is 

(t - T)^ = 2 . 45 X 10 ^ sin 


(91) 


6.6 TEHMS WITH PERIODS EQUAL TO STTs^ODiC PERIODS FOR JUPITER 
AND SATURN 


The eighth term, of Equation (44) is 



Terms will be obtained for Jupiter and Saturn due to their contribution to r„. 

. S . . . 

The contribution to r_ due to planet i with gravitational constant p. is 




(93) 


where pg is the gravitational constant of the sun augmented by the gravitational 

constants of the planets and moon. It is assumed that the helLccentric orbits 

of the earth-moon barycenter and planet i are circular and that the Latter Lies 

» S B 

in the ecliptic plane. With these assumptions, the vectors r^ and r^ , with 
rectangular components referred to the mean equinox and ecliptic of date, are 
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.s 

r . 

— 1 


s . 

1 


(94) 


r - sin, 

cos 'Ll. 
1 

0 



■ cos L ■ 
sin L 


a 



(9-5) 

i 


where s. is the circular orhit velocity of planet i relative to the sun (comjputed 
from Equation 89 with subscript J replaced, by i), and is the heliocentric 
mean longitude of planet i, referred to the mean equinox and ecliptic of date. 
Substituting (93) to (9'5) into (92) gives the contrihutioii to (t- t)g due to planet i: 


|i, s .a 

^ ~T^ (L-L^) (96) 

^ l^S 

Substituting numerical values from Appendix A for i = Jupiter gives 

(t-T)g^ = 20. 73 X 10’^ sin (L-Ej) (97) 


For i = Saturn, 


(t- 


'§SA 


4. 58 X 10 ^ sin (L-L 


SA' 


(98) 


6.7 DAILY TEIIMS DUE TO MOTION OF SUN RELATIVE TO SOLAR 
SYSTEM BARYCENTER 

The ninth term of Equation (44) is 



(99) 
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. • s 

Terms will be obtained for Jupiter and Saturn due to their eontribution to 
The eontribution to (t - r due to a planet i is obtained by substituting (93), 
(94), and (64) into (99) and using trigonometrie identities; 

[ 1 . s. 

(t - = - 2' 

- (1-cGs e ) sin (0^+\+L^) 

The second and third terms of this expression are neglected. Substituting (69) 
into the first term of (100) gives 

fi. s . ( 1+eos e ) 

(t- r)g^ = — - - ^ ^ ^ sin (UT1 + \^L-L.p (101) 


+ V sin e cos L. 


u 


(1 + e o s e) 


sin (0^ .+ \ 
M 


■L.) 


Substituting numerical values from Appendix A for i Jupiter gives 


{t-f)gj = 1-33 X 10"^^ u sin (UTl+X + L-Lj) 


( 102 ) 


For i ^ Saturn, 


(t- T)^g^ = 2.9 X 10"^^ u sin (UTl + X+L-Lg^) 


(103) 
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7, Final Exj)ressiQJi and. Arginine Jits 


Tike periodic terms of t - t, with analytLeal expressions for the 
coefficients, are given by Eqiaations (51 ), (58>, (71), (83), (88) (which applies 
for Jiapiter) and the corresponding term for Saturn, and (96) and (101), each 
evaluated for LTupiter and for Saturn. Substituting these terms, except the 
term of (83) with argument UTl + \ -f ZM, into (44) gives 


t - r = ATa t ^ 

e 


e sin E -t 


s a, . 
e M 

e^(l-tp) 


sin D 


i (l+eose)u f 


Ze‘ 




e sin (UTI+\) 


+ e sin (UT1+ X - M) +-| sin (UT1+ X - 2M) 


S (1 -GOS € )u 

[sin (UT14 \ + ZE) 

Ze"^ 


+ e sin (UTU- X + 2E -t- M;)] 


2e^ (l + ^i) 


u sin 


s (sin e }v 

+ X - p) - -2 eos L 


(104) 


® J . P , '"'SA ®SA . 

+ sm Ej 4 2 . SI® Eg^ 

^SA 


+ 


4 


sin(E-Lj) 4 ^ ^ sln(E-Lg^) 

G fJ-g 

S , (1+GOS € ) 

- — ^ — ... - u sin (UT14 X 4 L-L,) 

t^CA € ) 

- - z — U sin (UT14 X 4 L-Lg^} 

2g Pg' 


Substituting numerical values for the coefficients of the periodic terms, 
obtained from Equations (52), (60), (72), (84), (90), (91), (97), (98), (102), 
and (103), gives 
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t - T = AT^ + 1. 658 X 10'^ sin E + 1. 548 X 10“^ sin D 

+ 3.17679 X u sin(UTl+M + 5.312 X 16-12 ^ sin (UTl +\-M) 

+ 1,00 X lO'l^^g^j^ (UT19-K-2M} - 1.3677 X IQ-H u sin (UT1^X+2L) 

- 2.29 X 10“13 u sin (UT1+X+2E+M) 

+ 1.33 X 10-13 u sin (UTl-fK-D) - 1.3184 XiO-l^v gos L {105) 

+ 5. 21 X 10-^ sin Ej + 2. 45 X lO'^ sin Eg^ 

+ 20.73 X 10"^ sin (E-Lj) -5- 4. 58 X 10“^ sin (L-Eg^) 

+ 1.33 X 10-13 usin (UT14\4E-Lj.) 

•i 2.9 X 10"!^ usin (DT14X+E-Lg^) 

wkere tke e© ordinates n and v of tke atomiG clock wMck reads International 

Atomic Time t are in kilometers and t - r is in seconds. Tke coefficients of 

tke daily terms of (105) are proportional to n. For a clock at tke equator 

wkere n = 6378 km, tke magnitnde of tke coefficients of tke daily terms 

varies from a maximnm of 2. 026 X 10 s to a minimum of 1, 85 X 10 ^^s. 

Tke annual term witk argument E kas a coefficient proportional to v; its max- 

-6 

imum value is 0. 838 x 10 s. 


Tke last four terms of (105) are due to tke offset of tke solar 
system kar ye enter from tke sun. If tkese terms are deleted in (105)., coordi- 
nate time t in tke solar system karyeentrie frame of reference is replaced in 
tkis equation ky coordinate time t' in tke keliocentrie frame of reference, and 


tke acGuracy of tke equation is reduced sligMly. Tke previously derived ex- 
pression for t' - tis Equation (65) of Moyer (1971). In tkis previous expression, 
tke long -period potential terms due to Jupiter and Saturn ^tke terms of 105 
witk arguments E^ and were omitted along witk tke term witk argument 

E, whiek kas a significant effect on computed tkree-way range and doppler 
ofeservables. Tke previous expression for t' ~ t consisted of tke first nine 


terms of (105) witk minor differences in tke coefficients of tke periodic terms. 
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From Smart (I960), the relation "between the eecentrie and mean 
anomalies is: 


E = M -t e sin M + sin 2M + , . . (106 ) 

The three eGcentrlG anomalies in (ID'S) can "be Gompnted to snffieient acenracy 
from the following approximations to (106); 


E 

- M 4 e sin M 

(107) 


= Mj 

(10 S) 

^SA 

“ ^SA 

(109) 


The errors dne to the approximations in (107) - (109) are inGlnded in the error 
summary in Section §. Expressions are given "below for the arguments E, M, 


D, (L-Lj), (L-Lg^). 




which are fuii'Ctions of coordinate time t. 


The angle UTl in radians is eomipnted from universal time UTl in seconds 
using Equation (70), 


The Explanatory Supplement to the Ephemeris (1961) gives poly- 
nomials for the mean orhital elements E, M, and D (pp. 9S and 107). The 
argument for these expressions (quadratics or eu’bies) is Julian centuries of 
36 S2S days of ephemeris tim.e ET (coordinate time t) from January 0, 1900, 
12^ ET. Linear expressions have "been obtained which are tangent to the 
polynomials at January 1, 197S, 6^ ET (0.75 Julian centuries past the above 
epoch), hti units of radians, the linear expressions are 

L = 4. SS§ 339 1 1. 991 063 §3 X 10"'^ t 
M = 6. 24§ 291 + 1. 990 96S 71 x lo"”^ t 
D = 2.518 411 42.462 600 S18 X lo"^ t 

ItL 

where t is seconds of ephemieris time past January 1, 1950, 0^ ET, 

Seidelmann et al (1974) give quadratic expressions for the helio- 
centric mean longitudes of the earth, Jupiter, and Saturn (L„, L,, and L„. 

hj J 

respectively) and for the longitudes of perihelion for Jupiter and Saturn 


( 110 ) 

( 111 ) 

( 112 ) 
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(Wj and respectively), referred to tke mean eqmnox and eeliptie of 

1950.0. Linear expressions were olatained whiek are tangent to tke quadratic 

tl 

expressions for these mean orbital elements at January 1, 1975, 6 ET. The 

argument was changed from tropical centuries of 36 5E4. 2 19 §8 ephemeris 

h. 

days past 1950. 0 to ephemeris seconds past January 1, 1950, 0 ET. The 


expressions for the arguments L-Lj, L-Lg^, 


Mj, and Mg^ were obtained from 
the following combinations of the linear ej^ressions for L^, Lj, Lg^, j, 
and oj 


SA* 


L-Lj = 

^E- 

Lj ± ^ 

(113) 

I' - Sa = 

^E ■ 

Sa*" 

(114) 

Mj = 


O) 

J 

(115) 

•^SA = 

^SA 

" ^S.A 

(116) 


The resulting expressions, in radians, are 

L - L^ = 5.652 593 -t- 1. B23 136 37 X lo"”^ t (117) 

L - Lg^ s Z, 125 474 4 1. 923 399 23 x lo"”^ t (118) 

Mj - 5. 286 877 4 1.678 506 3 x 10 ‘® t (119) 

MgA = 1. 165 341 4 0.675 856 8 x lo"^ t (120) 


The me an longitude s on the left-hand side of (113) and (114) are referred to 
the mean equinox and ecliptic of date. The mean longitudes on the right-hctn4 
side of these equations are those of Seidelmann et al which are referred to 
the mean equinox and ecliptic of 19'50.0. The equinox difference does not 
affect the computed longitude differences L-Lj and L-Lg^, The difference in 
the reference planes does produee errors in Equations (117) and (118); how- 
ever, because of the small inclinations of the orbit planes of the earth, Jupiter, 
and Saturn te the ecliptic of 1950, 0, the errors are negligible. 
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Tke errors in. eomptited range and doppler ofeservaMes diie to the 
errors in. the linear approximations (110) - (112) and (117) - (120) are 
Gompletely negligiMe for the time period 1^50 - 2000. 
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8 . Effect o£ Retained and Neglected Terms of t - T 

The primary purpose of this section is to give estimates for errors in 
computed range and doppler observables due to errors in Equation (105), the 
final expression for t - T. The secondary purpose is to give the effect of the 
retained terms of t - r on these observables. 

Periodic terms of t T have a direct and an indirect effect on computed 
range observables. These two effects can be understood by considering the 
general procedure for comiputing a range observable. The light tirne solution 
produces solar system barycentrie position vectors of the transmitting station 
on earth at the transmission time t^ the spacecraft at the reflection time t^, 
and the receiving station on earth at the reception time t^- The reception time 
is known in atomic time r; addition of t - r from (105) converts it to coordinate 
time t (in the solar system haryeentrie space-time frame of reference). The 
light time solutidn yields the epochs and t^ in coordinate time. The epoch 
t^ in coordinate time is directly affected by the magnitude of the periodic 
terms of t - t evaluated at t^; the epochs t 2 and t^ in coordinate time are affected 
by the same amount to an accuracy of 3 or 4 digits (depending upon the rauge 
rate of the spacecraft). 

The values of t^ and in the XJTl time scale are also obtained. The 
former is not affected by the terms of t - r; the latter is affected but not sig-^ 
nifieantly. The indirect effect of a periodic term of t - Ton a computed range 
observable is due to the effect of the term on the epochs t^, t^, and t^ and 
hence on the position vectors computed at these epochs. The direct effect is 
due to the values of the term at t^ and t^ which appear explicitly in the equa- 
tioii for computing range observatles (Appendix B, Equation B4). Doppler 
observables are computed from the difference of two computed range observ- 
ables divided by the difference in the two reception times. The direct and 
indh-eet effects of a term of t - T on these two range observables produce a 
corresponding direct and indirect effect on the doppler observable. 

Appendi'-: B gives the formulation for computing two-way or three-way 
raoge and doppler observables. An analysis of these equations produced the 
approximate equations given below for the maximum direct effect of a term of 
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t - T on range and doppler observables. The direct effect of a term of t - T 
depends upon the form of the term. All of the retained terms and almost all 
of the neglected terms have one of the following forms. The terms can be 
divided into daily terms (period P 1 day) and long -period terms (P ranging 
from roughly half a month to several years). The daily terms contain the 
longitude of the atomic clock in the argument; the long-period terms do not. 
The coefficients of the daily terms are proportional to the spin-axis distance 
u of the atomic clock. The coefficients of some of the long-period terms are 
proportional to the distance v of tb . atomic clock north of the earth's equa- 
torial plane. The only neglected terms of t T which do not fit into these 
categories are some of the terms due to nutation, polar motion, and solid- 
earth tides. 

4. 

Liet6P and SP refer to the maximum possible value of the direct effect 
of a term of t _ t on range and doppler observables, ex,pressed as the traeking 
station to spacecraft range and range -rate, respectively. The effects of a 
daily term of t - t on two-way or three-way data are: 


6 P 

= Me 

(121) 

6 P 

II 

n 

(122) 


where M is the magnitude of the coefficient of the term and c is the speed of 
light. For three -way data, the spacecraft range p at which these maximum 
effects occur depends upon the separation in longitude of the transmitting and 
receiving stations; the separation can be adjusted so that the maximum occurs 
at any range. For two-way data, the maximum effects at a given range are 
given by Equations (121) and (122) multiplied by sin • Hence, the 

maximum two-way effects increase with range, reaching a peak at a range of 
43.2 astronomical units. 

The following equations give the effects of a long -period term of t - T on 
two-way data. They also give the effect of a long-period term on three-way 
data if the coefficient of the term is not proportional to v. 


4 4 
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(123) 

(124) 

If the eoeffieieat of a long -period term is proportional to v, the effects 
of the term on three-way data are given by (121) and (122), where M is the 
maximinm possible magnitude of the coefficient of the term. These maximum 
effects are independent of the range to the spacecraft. 

The indirect effects of a term of t - T on two-way or three-way range and 
doppler observables will not exceed the following approximate values: 

- pM (125) 



SPj. = aM (126) 

where p is the tracking station to spacecraft range rate, and a is the Inertial 

acceleration of the spacecraft. The highest spacecraft acceleration likely to 

2 

be eneountered is 25 m/s , which occurs near the surface of Jupiter. Higher 
accelerations occur closer than 3.3 solar radii from the center of the sun, 
hut it is unlikely that a spacecraft would enter this region. 

An extensive error analysis has heen performed to determine estimates 
for the maximum possible errors in computed range and doppler observables 
due to errors in Equation (104) or (105) (which contain analytical and numeri- 
cal coefficients, respectively) for t - r. The results of this analysis are pre- 
sented in Table II. The second coilumn lists the errors committed in the 
derivation of Equation (105), If an error is associated with a term of (104) 
or (105), that term is identified in column 1. The magnitudes of the eoeff ieients 
of the error terms are given in column 3. The maximum values of the direct 
effects of the error terms on computed range and doppler observahles are 
given in columns 4 and 5, respectively. When one figure is given for the 
effect of an error, it applies for two-way or three-way data. When two 
figures are given, the figure in parentheses is the error for three-way data 
and the other figure applies for two-way data. For daily terms, the maximum 
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TABLE II 
Error sumim ary 


T e rm 

Sonar ce of error 

Magnitude of 
error term, 
p-s 

Direct effect on 
range, 
m 

Direct effect on 
range rate, 
pm / s 

Arguimie nl 
E 

Assmmlmg e cO'nstan't 
from 195 0 to 2 0i0'0 

1. 0 

1, 5 

0. 3 

Arguime 

E 

2 

N'e.g Le G te'd e te rm 
in Ei^uatioin {107) 

0. 3 

0. 7 

0. 3 

Argument 

E 

Ignorinig periodic variatioms 
in heLiocentric orbiLtal 
eLe meats of earth -moon 
ha ry center 

12. 0 

15. 0 

4. 0 

Argume nt 
D 

Ignorinig and e 

0. 2 

2. 3 

10, 0 

CoefflGients 
with factor 
s ^ (1 + e os e } 

Ignoring terrois 

0. 0'0 

0. 01 

0. 8 

Coefficients 

with factor 

s { 1 - G OiS e ) 
e 

2 

Ignoring e terms 

0. 00 

0. 02 

1. 2 

ALL with 
UTl + X in 

Appr o xi ma ti on i n 
Eqinatiom (69) 

0 . O'0 

0. 0® 

5. 5 


ar g lame nl 
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TABLE II 


Error siurrunary (eomtd) 


Term 

S'O'nrce of error 

Magnitude of 
error term, 
ps 

Direct effect on 
range, 
m 

Direct effe 
r a n ge r a 
pm /s 

Argume nt 
UTi + \ f 2M 

Deletimg this berm fr’om 
Equ-atioin (84) 

0. 00 

© . 02 

1 . 5 

All with 
factor s^iit 
or S(>v Ln 
coefficient 

Ignorinig periodic variatiors 
in he lioeentrie orhital 
eLemeats of earth-m'OOin 
baryee Piter 

0. 00 

0. 4 
(0.6) 

30. 0 

Ar'gnrnent 

UTI + \ - D 

Ignoring ej^ and 
second term of 
Equation (45) 

0. 00 

0. 03 

1. 9 

ArgiH'mie nt 

Ignorinig third term 

0. 00 

0 . 01 

0. 02 

UTI f X - D 

in Eq^aatioin (65) 


(0. 11) 

(0. 3) 

Argvvmie nt 

Ignorinig e 

0. 02 

0. 04 
(4.2) 

0, 02 
(1.7) 

Argnmients 
Ej a„,d 

ISfe glee ted e terms in 
E q^na 1 1 on s ( i 0'8 ) and' ( 1 09 ) 

0 . 2 

0. 04 

0. 010 


N-e g ie e te d te rm s i n 
Eq nations (37) and (39) 

1. 3 

1. 1 

0. 3 


Neglected gravltatioinal 
potential at earth doe to 
moon in Equation (27) 

0,. 00 

0 . 1 

0. 2 


■~a 
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TABLE II 

Error summary {eontd} 


Term 


Magnitude of Direct effect oa Direct effect on 

error term, range, range rate, 

Source of error jj.s m (j.m/s 


Arguments Igmorinig e, e_, and 2. 1 2,3 0. 5 

L - Lj and 
L - LgA 


Argument 


L - L. 

1 


Argument 
UTl + ^ + L - 
L 

1 


Ignoring terms for 
M'e r e ur y , V e nua , Mars, 
U r anus., Ne ptune , an d 
Pluto, including effects 
due to eecentrieities 


Omitting 

i /r^. 

c2 V — S 


the 

B\ 

^ e ) 


te rm 
in 


Equation (40) 

Ignoring secomd and third 
terms of Equation (1O|0) for 
all planets, the first term 
for all pilanets except 
Jupiter and Saturn, and 
ignoring eccentricities of 
planetary orbits for retained 
and ne glee ted terms 


1. 3 


0. 0'O 


0. O'0 


1.8 


0 . 02 


0 . 06 
{0.2) 


0. 4 


0. 05 


3. 6 
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TABLE II 


Error siamroiary (Gontd) 


Term 

So'iarce of error 

Magnitude of 
error te r m , 

ps 

Direct effect on 
range , 
m 

Direct effect on 
range rate, 
pm /s 

Argurriie Pits 
B, UTl + K- B, 

L - Lj, L - Lg^, 
UTl -1- X + L - 
Bj, anid UTl -f-\ 

+ L - LgA 

Ignoring ine Linatiora of 
orbits of mioon, Juipiter, 
and Satnrn to ecliptic plane 

0. 02 

0. 2 

1. 5 

Coefficients 
with factor 
1 -f- cost, 

1 - eos £ 
or sine 

Assuming t and e are 
constant from 195 0 to 
ZO'OO 

0. 00 

0. 02 
(0. 06) 

1.4 


Ignoring nutation, polar 
rmo ti on , a nd S'0 1 i d - 
earth tides 

0. 00 

0. 1 

3. 0 

All 

Roundoff of coefficients 

0. 2 

0. 3 

0. 7 


Misce llaneous 
{see text) 

4. 3 

3. 8 

7. 8 


Totals 

Z3. 0 

30. 0 
(34.6) 

75. 0 
(77. 0) 


v0 



values of the errors are listed. This is also true for three-way errors due 
to long-period terms which have coefficients proportional to v. The errors 
in two-way data due to long-period terms and in three-way data due to long- 
period terms whose coefficients are not proportional to v were esmputed for 
a spacecraft range of 50 astronomical units. 

A large numfeer of the errors listed, in Table II are due to terms 
neglected in (104) and (105). Another large group of errors is due to 
ignoring eccentricities of elliptical orbits or ignoring terms above a certain 
power of the eecentrieity. Some of the largest errors are due to ignoring 
periodic variations in the heliocentric orbital elements of the earth -moon 
baryeenter. These variations, in the form of longitudinal and radial per- 
turbations, were obtained from Newcomb (1898). The effects of periodic 
variations in the geocentric orbital elements of the moon and in the helio- 
centric orbital elements of Jupiter and Saturn have not been analyzed. The 
miscellaneous error listed at the end of Table II is an allowance for these 
errors and other minor unanalyzed errors. 

Estimates have been obtained for the maximum possibile direct errors 
in computed range and doppler observables as a function of the range to the 
spacecraft. For two-way data, these upper limits are 

68^= {0. 62 m) AU (127) 

= (2.4 X 10'^m/s)AU (128) 


where AU is the range to the spacecraft in astronomical units. These figures 
are based upon the constant derivatives of long-period errors with respect to 
range and the derivatives of the daily errors with respect to range evaluated 
at a range of zero. Since daily errors are prooortional to the sine of a 
multiple of the spacecraft range, these upper limits are high at large space- 
craft ranges. For three-way data, the upper limits are 

dPg = (0.59 m) AU + 5. 5 m (12 9) 
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(130) 


= (0.49 X 10'^ m/s) AU + 53 X 10"^ m/s 

From Table II, the estimate of the maximum error in the magnitudje of 

t - T is 23|±s. Using this figure and a spacecraft range rate of 30 km/s, 

Equation (125) gives a maximum indirect error in computed range observables 

2 

of 0.69 m. For a spacecraft inertial acceleration of 25 m/s , (126) gives an 

upper limit to the indirect error in computed doppler observsibles of 575 X 

“ 6 

10 m/s. This error is almost eight times as large as the maximum value 
of the direct error. 

To an accuracy of 0. 01 ps, the error in the magnitude of t - r is due 
entirely to long-period terms. This error could be represented by a curve 
which is a function of a small number of parameters. When fitting to tracking 
data obtained from a spacecraft which is near a planet, the values of these 
parameters could be estimated, thereby eliminating most of the indirect 
errors in computed observables due to neglected terms of t - T . 

Table III gives the same information for the retained terms of t-T as given 
in Table II for the neglected terms. It is seen that the maximum value of the 
direct error in a computed range observable due to neglected terms of fc - t is 
about 1% ©f the maximum value of the direct effect of the retained terms of 
t-T. The corresponding figure for a doppler observable is 0, 16%. The maxi- 
mum error in the magnitude ©ft - T is about 1. 36% of the maximum value of 
the sum of the periodic terms of t - T. Hence, the maximum value of the 
indirect error in a computed observable is about 1.36% of the maximum value 
of the indirect effect of the retained terms. 

Approximate maximum values for the direct effects of the retained 
terms of t - r, as a function of the spacecraft range AU (in astronomical units). 


= (74.3m)AU (131) 

= (1. 72 X LO^^m/s) AU (132) 
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TAILE III 


Maximum effect of retained terms 


Argument 
of term 

Magnitude 
of term 

HLS 

Direct effe ct 
on range, 
m 

Direct effect 
on range rate, 
[j.m / s 

E 

1658. 

2476. 

493. 

D 

1. 548 

28. 6 

70. 4 

UTl -1- X 

2 . 026 

607. 4 

44173. 

UTl + X - M 

0. 0339 

10. 2 

739. 

UTl + X-. 2M 

0. OO'O 64 

0. 19 

13. 9 

UTl + X + 2L 

0. 0872 

26. 2 

1902. 

UTl -f X -f 2L + M 

0, 0015 

0. 44 

31. 8 

UTl +X - D 

0 . OO'O 8 5 , 

0. 26 

18, 5 

L 

0. 838 

1. 25 

0. 25 



{251. ) 

(50. 0) 


S. 21 

0. 66 

0. 01 

^SA 

2. 45 

0.13 

0. O'O 


20. 73 

28, 4 

5. 17 

^-^SA 

4. 58 

6. 61 

1. 27 

UTl + X + L - ,L 

u 

0. OO'O 85 

0. 26 

18. 5 

UTl + X+ L - Ug^ 

0. OO'O 19 

0. 06 

4. 03 

T otals 

1696. 

3187. 

47 471, 



(3437. ) 

(47 521. ) 


2 


JPL, Teekmeal Memoramdium 3 3“ 7 86 


(51ed) AU + §% cn 


(133) 


6 P, 


aPg = (0. 012 X lQ"^m/s) AU + 47. 0 X 10"^ tn/s 


(134) 


2 

For a sg)aeeeraft range rate of 30 km/s and inertial acceleration of 25 m/s , 
the maximu;m valiaes of the indirect effects of the retained terms of t - r on 

. -5 

computed range and doppier observables are 51 m and 42 x 10“ m/s, 
respective Ly . 
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9. ALfeernate Exf>yessioJi for t T 


This section gives an alternate eqnatioh, for fc- T which is a fanction of 
position and velocity vectors of the atomic clock and the major bodies of the 
solar system. The eqiuation is obtained from Equation (44) by replacing 
the second, sixth, and seventh terms, which are expressed as integrals, 
with fnnetions of position and velocity vectors derived in this section. 

The second term of Eqination (44) is given by Equation (SI). Eor the 
heliocentric elliptical orbit of the earth -moon barycenter, 

, S S j ' 

“ *^1 ' ~ a e sin E 

Snbstitnting the approximate form of (135) into (51) gives the desired form 
for the second term of Equation (44); 


(t - r)2 . 


2 / .S S\ 

2 • *^b) 


( 136 ; 


The sixth term of Equation (44) is given by Equation (§§). The equation 
analogous to (135) which applies for the heliocentric orbit of Jupiter is; 




(137) 


Substituting e ^ sin Ej from (137) into (88) and using (89) gives 


(t - t), 




,.S S, 


The corresponding term for Saturn (SA) is the seventh term of Equation (44); 
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( 139 ) 


JPL 


( t - t) 


- gA ,.S S , 

^ 2 . (£g^ • £g^) 




Replaeidg terms two, six, and 


even of Eqnntion (44) with Equations 


(136). (138) aad (139). the three terms which contain rf. and 

reordering the terms gives the alternate expression for t - r: 


t - r = at 4 . i t-C El 1 .s 


E ■ -A^ 

G 


1 t-C s 

% • 




/•S s, 

^ 2 


*^SA ...^ 


(E-ta . £ 


SA' 


SA ^SA' 


If the term (rS' . r^) 
_2 '-S 


.S . 


(140) 


, which was omitted in Equation (40), 


was 


reinstated, f ^ in the fourth term of (140) would be replaced hy f‘=. 
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Appendix A: Notation and N^rneriGal 


Values for Farameters 


This appendix gives the definitions for the parameters used globally 
throughout the text. Numerical values are given for those parameters which 
appear in the final expression for t ~ t. 


Subscripts and Superscripts 

A = location of atomic clock on earth w" ■ 'h reads International 
Atomic Time t 


E = earth 

B = earth -moon bar yc enter 

M = moon 

S = sun 


C = solar system baryeenter 

J = Jupiter 

SA = Saturn 

p This subseript indicates that periodic terms of the quantity are 

to be retained and constant terms are to be discarded 

o This subscript indicates the quantity is evaluated at the initial 


A bar over a quantity indicates the time average value of the quantity. 

Position, Velocity, and Acceleration 

J"!* = position, velocity, and acceleration vectors of point i 

relative to point j. The dots denote differentiatipn with 
respect to coordinate time t 

s^ = velocity of point i relative to point j 

I = velocity of fixed atomic clock on earth relative to solar system, 

baryeenter 


5b 


SEPEODUCBTLITY OF Til* 
PAGE IS FOOE 


JPLi Technical Memorandum 33-786 


Time 


T 

T-i- 


t 


n 


= proper time obtained from an atomic clock on earth 

= International Atomic Time (TAI) obtained from an atomic clock 
on earth. Starting in Section 4 , T* is denoted by t 

= coordinate time in solar system baryeentric space -time frame 
of reference 

= conversion factor from cycles obtained from an atomic clock 
to seconds of atomic time t 


n-!' = conversion factor from cycles obtained from a cesium atomic 

clock to seconds of International Atomic Time 

AT^ = constant term in expression for t - 


UT 1 = observed nniver sal time, corrected for polar motion. Equa- 

tion (70) converts UT 1 from, seconds to radians. 


Physical Constants 


The physical constants were obtained from. Standish et al (197b) or were 
computed from quantities obtained from this reference. 


e 


A 


E 


= speed of Light - 29'9 792.458 km/s 

= the number of kilometers per astronomical unit (AU) 


= 149 597 871.410 5b km/AU 

3 2 

= gravitational constant of body i, km /s 

= Gm^, where G is the universal constant of gravitation and m^ 
is the mass of body i. 


The gravitational constant of the sun is computed from A^ using Equation (104) 
of Moyer (1971): 




S 


E 


(86400)^ 


1. 32712442 xlG^^ km^/s^ 
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where 


k 


= the Gaiissiati gravitational eonstant 
= 0.Q17 202 098 95 AU^^^/day (exactly) 


Using ^ and. mass ratios ofetained from Sfeandish et al (1976) gives 
o 

fjL^ = 398 60'O. 6 km^/s^ 

= 4 902.79 km^/s^ 
laj = 1.267 120X10^ km^/s^ 

*^SA ^ 3. 793 410X10'^ km^/s^ 


|4 



§1 ,3 00 7 




= gravitational eonstant of the sna augmented by the gravitational 
Gonstants of the planets and the moon 

= 1 . 32S 906x10^^ km^/s^ 


Gravitational Potential 


U = Newtonian gravitational potential, using the positive sign con- 
vention (i.e., U = -<f>) 

U. = U at location i 
1 

U(j) = U due to body j 

U^(j) = U at location i due to body j 


(VU)^ = gradient of U at location i 


Station Location 

The earth -fixed coordinates u, v, and X. of the atomic clock which reads 
International Atomic Time are: 


u - distance from earth's spin axis, km 


JFL Technical Memorandum 33-786 



V 


- distance north of earth's equatorial plane, ktn 
X - east longitude 

®M " cnean sidereal time = Greenwich hour angle of mean equinox of 
date 


Meliocentrie Or hit of Earth -Moon Barye enter 

Numerical values were obtained from page 98 of the Ex.planatory Sup- 
plement to the Ephemeris (1961). 


= semi -major axis = 1.000 000 23 AU 
= 149 §97 906 km (obtained using A given above) 

ill 

= eceentrieity = 0.0167 2. From 1950 to 2000, the last 
changes from 3 to 1 . 

= radial coordinate 

= velocity 


- circular orbit velocity 


M 

E 

f 

L, & 


V 


e 


= 29. 784734 km/s 

= mean anomaly 

= eccentric anomaly 

= true anomaly 

= mean and true longitudes, respectively, of the sun, measured 
at the earth-moon barycenter. These angles are referred to 
the mean equinox and eeliptie of date. 

= elevation angle of velocity vector, measured from the transverse 
direction (normal to radius) 

= mean obliqiuity of the eeliptie 

- inclination of ecliptic plane to mean earth equator of date 
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cos e 


0. 917 46 (frocn 1950 to 2000, the last digit chatiges from 4 to 8) 


sia e = 0.397 83 (from 19'50 to 2000, the last two digits chaage from 
88 to 78) 

He lie centric Orbit of Planet i 

a^ = semi-major axis. From Seidelmann et al (1974), j: 

aj = 5.202 833 481 AU = 7. 783 32 8 x 10^ km 

^SA " 9. 538 762 055 AU = 1 . 426 978 x lo"^ km 

= eccentricity . From Seidelmann et al (1974), evalmated at 1975.0, 

ej = 0, 048 2 84 

e^A = 0.056 038 

These figures are constant, to three significant digits from 1950 to 

2 000 . 

r, = radial coordinate 

L 

s^ = e ire alar orbit velocity 



Sj = 13.064 13 km/s 
~ 9-645 16 km/s 
= mean anomaly 
“ eccentric anomaly 

Lj^, = heliocentric mean and tme longitiades, respectively, of 

planet i, referred to the mean equinox and ecliptic of date. 
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Geocentrie Orbit of Moon 

= semi-tnajor axis = 384 3 99. 1 km 

This value is obtained from the above values of and and the observed 
mean motion of the moon using Equation (106) of Moyer (1971), which is a 
modified version of Kepler's third law. 

= eccentricity = 0. 054 9 (used in error analysis only) 

= circular orbit velocity 

= 1 . 024 548 km /s 

J = geocentric mean longitude of the moon, referred ho the mean 
equinox and ecliptic of date 

D = C - L = mean elongation of the moon from the sun 
Miscellaneous 

(t - r)j. = term i of Equation (44) 
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Appendix B; CQoapia.tatlon. of Range and Doppler Observables 


This appendix gives a barief description of sicnplified versions of the 
formulas used to obtain the computed values of range and doppier observables 
obtained by the Deep Space Network of the Jet Propulsion Laboratory, The purpose 
is to show how the time transformation t - r is used in the computation of these 
observables. The simplifications made to the formulation do not change the 
effects of t - T. For further details of the formulation., see Moyer (1971). 

The eomiputation of two-way or three-way range observables Is des- 
cribed first. A signal is transmitted from a tracking station on earth at 
time tj^, received and retransmitted at the spacecraft at time t^, and received 
at the samie tracking station on earth (two-way data) or at a different station 
(three-way data) at time t^ . The definition of the range observable R is 

R = tg (TAI) - tj (TAI) (Bi) 

where t^ (TAI) and t^^ (TAI) are the receptio-n and transmission times, respec- 
tively, in Inter national Atomic Time TAI. The "time tag" associated with 
each range observable is the known reception time (TAI). 

The first step in eomputing a range observable is to obtain the light 
time solution. The epoch t^ (TAI) is converted to coordinate time in the solar 
system barycentrie space-time fram.e of reference (denoted here as ephemeris 
time ET) using 


t3(ET) - t3 (TAI) -t- (ET-TAI)j. 


(B2) 


me tignt 


where (ET-TAI) is (t-r) given by Equation (lOf), evaluated at t 
time solution gives the epochs (ET) and t, (ET) and the position vectors 

c c c ^ ^ 

(t^), £2 (^ 2 ^* —1 ^^1^’ solar system barycentrie position 

vectors of the receiving station at t^, the spacecraft at t^* and the transmit- 
ting station at t^, respectively. The epochs t^ and t^^ are also obtained in the 
UT 1 time scale; they are used along with the ET values of these epochs to 
compute the geocentric position vectors of the receiving and transmitting 
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stations at t^ and t^, respectively. These vectors are used in forming the 
above solar system barycentrie vectors. The following vector magnitudes 

are eomjputed: 


4= (ti) 


i = 1,2, and 3 


12 


^23 


42 -4i 

^3 (^ 3 ) - (t^) 


^2 '■'2 


The epochs t^, t 2 » and t^ satisfy the following equations; 


23 


■P-<s 


t^iET) - t2(ET) = 


Hn 


^2'*'^3’‘'^13_ 


^ 2 +^ 3"^23 


3^2 

2^1 


(B31 


The first term on the right-hand side is the Newtonian light time; the second 
term is an approximate expression for the contribution to the light time 
from general relativity. 


Given the Light time solution, the range observable is computed from 





G 



fn 


\"l+^2-^12/ 



3 


£n 


^ 2 +^ 3 -" 23 / 


- (ET-TAI) + (ET-TAI) + A 


(B4) 


The sum of the first four terms is the round-trip light time in ET (see 
Equation E3). The next two terms convert this interval to an interval of TAI. 
The last term includes corrections for the effects of the troposphere and 
charged particles. If the light time solution has a tolerance 6t for the deter- 
mination of t 2 and t^, observables computed from (B4) will have a corres- 
ponding error of roughly (tf ^2+^23^*®*''^^' ^ i*ange observables were computed 

from (Bl) using t^(TAI) obtained from the light time solution, the error 
would be about dt, which is approximately four orders of magnitude Larger. 
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It is seen that the titme ti’anforteiation (t - r) affeets range observables 
directly through terms 5 and 6 of (B4), it also affeets the epochs t^(ET), 
t 2 i(ET), and t^(ET) and hence the quantities ^ 23 * ^2' ^3 

this is the indirect effect of (t - t). 


Doppler observables are derived from a signal being continuously' 
transmitted from, the transmitting station on earth to the receiving station via 
the spacecraft. A particular observable is associated with an interval of 
reception (in seconds of TAI) at the receiving station. Typical values of 
are bO s and 600 s. It can be shown that the value of a dioppler observable 
obtained by the Deep Space Network is equal to 


F 


IT 

f 


(a -R ) 

' e s 


where f™ is the frequency of the transmitted signal, and R and R are pseudo 
round-trip range observables, defined by (Bl), with reception times equal to 
the end (e) and start (s), respectively, of the reception interval T^. Com- 
puted values of two-way or three-way doppier obaervabies are obtained from 
(B5) using pseudo two-way or three-way range observables eO'caputed from. 
(B4). The computed doppier observables are in error only because of errors 
in the computed pseudo range observables. 


When (B4) is used to obtain the computed value of a true range observ- 
able, the contribution to A due to charged particles is positive because the 
ranging signal travels at the group velocity (<c). When pseudo range observ- 
ables are computed (which are differenced to form doppier observables), the 
sign of the charged particle correction is negative because the doppier signal 
travels at the phase velocity (>e). These two eorreetions, which have opposite 
signs, have the sacne magnitude. 
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